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Abstract: 

Matter implies the existence of a large-scale connected cluster of a uniform nature. 
The appearance of such clusters as function of hadron density is specified by percolation 
theory. We can therefore formulate the freeze-out of interacting hadronic matter in terms 
of the percolation of hadronic clusters. The resulting freeze-out condition as function of 
temperature and baryo-chemical potential interpolates between resonance gas behaviour 
at low baryon density and repulsive nucleonic matter at low temperature, and it agrees 
well with data. 



Consider a hot quark-gluon plasma in thermodynamic equilibrium, specified in terms of 
temperature T and baryochemical potential /i. Reducing the temperature of this medium 
at constant fi eventually brings it to the confinement transition at Tc,^{Tc). Below this 
temperature, the system consists of interacting hadrons. Further cooling finally leads to 
freeze-out, at Tj, fi(Tf); beyond this point, we have non-interacting hadrons. The resulting 
two dividing lines in the T — fi plane are schematically shown in Fig. ^ 
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Figure 1: States of matter in QCD - general view 
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While the hne separating the deconfined QGP state from the confined hadronic medium 
is calculable in lattice QCD at finite temperature and density (although so far with consid- 
erable problems at large fi and low T), the freeze-out curve is less well-defined theoretically 
as well as experimentally. The aim of the present paper is to address how freeze-out can 
be specified conceptually. 

In the two limiting cases of hot hadronic matter of vanishing overall baryon density and 
cold nuclear matter at vanishing temperature, we can make use of some specific aspects 
of the relevant dynamics. The confinement transition leads to an interacting hadronic 
medium. Now it is known that if the interactions between the constituents of such a 
medium are dominated by resonance formation and decay, the system can be treated as 
an ideal gas of all possible resonances ^1 12] • In our context this means that if and when 
the interacting hadronic medium is resonance dominated, freeze-out effectively occurs at 
the point of confinement. In particular, the relative abundances of the different species of 
hadrons and hadronic resonances are in this case determined at the transition from QGP 
to hadronic matter. 

For systems of vanishing or low baryon density, this appears quite well supported. On 
a theoretical level, it is in fact claimed that resonance formation dominates the interaction 
between hadrons, perhaps most clearly in the dual resonance model [3]. Experimentally, 
an analysis of the species abundances indicate that these are indeed determined by a 
single freeze-out temperature, Tf ~ 175 MeV, obtained e~^e~ annihilation, p — p and p — p 
interactions as well as in heavy ion collisions P| . This temperature is moreover completely 
in accord with that found for the confinement transition in finite temperature lattice QCD 



On the other extreme, for dense nuclear matter at low temperature, the situation is 
quite different. The interaction between two nucleons does not lead to resonance forma- 
tion; instead, it is dominated by Fermi statistics and baryon repulsion. In particular, 
at T = 0, freeze-out must occur at a density equal to that of normal nuclear matter, 
riQ ~ 0.17 fm~^. If we consider a system of non-interacting nucleons at T = 0, the 
corresponding value of the baryochemical potential is obtained from 

no(T = 0, /i) = -^(/.^ - m2)3/2 = o.l7 fm^ (1) 

and found to be /iq — 0.979 GeV. Note that for T = and this value of fif, the system 
consists only of nucleons. The threshold for the occurrence of baryon resonances is /i = 
1.23 GeV, corresponding to 4.2 fm~^ or about 25 times normal nuclear density. Hence 
interacting cold nuclear matter very likely contains only nucleons over the entire range 
from freeze-out at /x/ = /io GeV to the deconfinement value /ic, where it becomes cold 
quark matter. The freeze-out value of fi thus obtained is, however, in principle somewhat 
model-dependent. We have here assumed the repulsion to be totally given by Fermi 
statistics; there could well be additional repulsion beyond this. Using hard core repulsion 
with a core volume Vq |6^ , Eq. [T] is replaced by 

frj. n A no(T = 0,/i) ni7f 3 fo\ 

nviT = 0, /i = — — — = 0.17 fm . 2 

1 + ?2o(T = 0,/i)Vo 

The result is a freeze-out at slightly higher values of yU. To illustrate, for Vq = (47r/3)-R^ 
and a nucleon radius of 0.8 fm, freeze-out occurs at [ly — 0.992 GeV. Such a hard core is 
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presumably too large; reducing R will reduce the resulting /iy, so that the physical value 
falls most likely between the values of /io and fiy 

Empirically, the determination of freeze-out parameters through the relative abun- 
dance of hadron species so far appears to be the only unambiguous approach to the 
problem. We therefore adopt this "chemical" freeze-out definition in the remainder of 
the paper. As an immediate consequence, we have to revise Fig. [T] At /z = 0, T/ = Tc, 
while at T = 0, /i/ ~ /io < /Uc, as shown in Fig. |21 From what was said above, it is clear 
that for very low temperatures, specifying freeze-out through species abundances becomes 
problematic, since near T = 0, the system contains only nucleons. This has to be kept in 
mind when studying cool baryon-rich media. 
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Figure 2: States of matter in QCD, including resonance gas effects 

Looking at Fig. |21 we note that we have so far identified only Tj = Tc at fi = and fij 
at T = 0. To determine the entire freeze-out curve in the T — fi plane, we have to know 
how the contributions from non-resonant baryon interactions modify the ideal resonance 
gas picture, in order to specify how the freeze-out curve departs from the confinement 
curve. In other words, we have to identify a freeze-out condition. 

A very successful phenomenological parametrization of the available heavy ion 
data is obtained by requiring the average hadron energy per average hadron number to 
be constant at freeze-out, 

(E) 



(AT) 



1 GeV. 



(3) 



While this parametrization accounts well for all data from SIS to RHIC experiments, it 
is clear that it encounters difficulties for sufficiently low temperature. At T = 0, it leads 
to freeze-out at /x ^ 1.039 GeV. With Fermi repulsion only, this corresponds to ~ 
0.73 fm~^ ~ 4.6 no, with the hard core repulsion used above it gives nj ~ 0.285 fm^^ ~ 
1.7 no. Both values would exclude the existence of normal nuclear matter, for which 
(E) / (N) is below 1 GeV. Since it is also not known what features of the underlying 
physics could lead to Eq. (fT^ as freeze-out condition, it seems worthwhile to study possible 
freeze-out mechanisms and attempt to find a consistent description for the entire T — /i 
plane. 
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Matter implies the existence of a large-scale interconnected medium of uniform nature. 
When such a system breaks up into fragments much smaller than the size of the volume 
in which it is contained, it has undergone a change of state. One possible way to define 
freeze-out is thus geometric: it occurs at the point at which the size of the largest hadronic 
clusters falls below the size of the given overall spatial volume. This point is determined 
in percolation theory jH] and for three space dimenensions becomes 

Vh 

where n = N/V specifies the hadron density, with hadrons in the overall volume V. 
The volume of an individual hadron is denoted by Vh = (47r/3)r^, and hadrons are allowed 
to overlap; thus Vh introduces the short range nature (~ Vh) of hadronic forces |9j. 

When the hadron density has reached the percolation point n = nf, the largest clusters 
reach the size of the overall volume. This, however, does not imply that they fill the 
volume. In fact, it is known that at the percolation point, still exp{— 0.34} ~ 71% remain 
empty space. Hence the vacuum, measured in terms of the hadronic scale Vh, also forms 
percolating clusters. We can therefore ask for what density vacuum percolation stops and 
only the strongly interacting medium spans the entire space. This occurs for 

nc = -r^, (5) 

obtained in an analogous way as Eq.^ For densities above n = ric, any large-scale vacuum 
has disappeared. Since the disappearence of the physical vacuum is a basic feature of 
deconfinement, it seems natural to relate this threshold to the confinement /deconfinement 
transition. 

Hence, on a purely geometric basis in terms of connected clusters, we have two thresh- 
olds: a.t n = Uf, there appear clusters of the size of the overall spatial volume, and at 
n = ric, the vacuum disappears as a connected medium 0. Let us see if and how they be 
related to the known freeze-out points at T = and at fi = 0. 

At T = 0, we expect freeze-out to occur only when the nucleons no longer form 
interconnected matter. If we use the nucleon radius r„ ~ 0.8 fm for Vh, Eq. (j3)) leads to 
Uf ^ 0.16 fm~'^ as the freeze-out density, and this is in good approximation the density of 
normal nuclear matter. The correspondence between hadron percolation and freeze-out 
thus works well at T = 0. It defines nuclear matter as the most dilute system of nucleons 
which still forms connected matter. 

Can we also relate the other extreme, freeze-out and deconfinement at /i = 0, to 
percolation? Eq. Q determines the density of hadrons for which the vacuum disappears 
as a connected medium. From the arguments given above, we can consider the interacting 
hadron system at /i = as an ideal resonance gas. Hence it seems natural to use the 
vacuum percolation condition 

1.24 

(4^rr 

for the ideal resonance gas density n{T,fi = 0) to determine both freeze-out and decon- 
finement. To obtain n{T,^) and other necessary observables of such a system, we briefly 
recall the essentials of the ideal resonance gas model. 



n{T,^l = 0) = JJ-J:^, (6) 
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In thermodynamics, a closed system of fixed energy E, volume V, and number Ni 
of particles of type i is described as micro canonical ensemble. Bringing this system in 
contact with a heat bath leads to the canonical ensemble, having an average energy 
(E) = tr(pif), while V and Ni are conserved exactly. If we further also let the particle 
numbers Ni fluctuate with only their averages conserved, we obtain the grand canonical 
ensemble. For this, we have 

where Aj = e^'^* is the fugacity of particle species i and Z^. the corresponding canonical 
partition function. The chemical potentials fii assure particle number conservation in 
an average sense. Changing the summation over discrete quantum states to phase space 
integration, we obtain the grand canonical partition in the form 



In Zg = Y1 9i7r^ / dpp'^ In 1 + 77^ AiC 
^ 2tt^ Jo l 



Vi 
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where p° = ^p'^ + mf, with rui for the mass of the particle species i, Qi counts the 
spin/isospin degeneracy, and rji is the statistics factor, with rji = 1 for fermions and 
rji = —1 for bosons. 

In a relativistic gas, in which the particle numbers are generally not conserved, the 
chemical potentials are associated to the conserved baryon, charge and strangeness quan- 
tum numbers B, Q and S. The abundance of strangeness is complicated by the larger 
strange quark mass, which could prevent that strange quarks are produced as abundantly 
as the lighter non-strange quarks. In order to allow for a resulting reduced strangeness 
production, the parameter 7^ (7^ < 1) was introduced [Hlj. For 7^ = 1, no suppression 
occurs, while for 7^ < 1, fewer strange quarks appear, with a corresponding suppression 
of strange particle production. A factor 7^ is associated to each strange quark, so that for 
hadrons containing fg strange quarks, there is an additional multiplier 7/* in the fugacity. 
We thus have 

A. = 7/^AfA«'Af , (9) 

where Bi, Qi and 5*, are the quantum numbers for individual particle species. Then the 
particle numbers are given by 

{N,)=g,— rfpp^ (7/^A|»A§>Af)-VP^+d" , (10) 



'27r2 Jo 

and the overall density of hadronic constituents by 

niT,^^) = ^p^. (11) 

The average net baryon number, i.e., the number of baryons minus that of antibaryons, 
is the sum over the particle numbers weighted by the baryon quantum number, 

{NB) = j:B.m: (12) 
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and the corresponding baryon density is 



nsiT, /i) 



(Nb) 
V ' 



From 



we obtain with 



V ' 



(13) 
(14) 
(15) 



a similar relation for the energy density. 

For /X = 0, we now use Eq. (fTT|) in the percolation condition (0), including all ob- 
served resonances up to mass 2.5 GeV, with the nucleon radius Vh = 0.8 fm for baryons; 
for mesons, we use a value smaller by a factor (2/3)^/^ ^ 1.1, as suggested by bag model 
arguments. This determines the threshold temperature for confinement and (species abun- 
dance) freeze-out at vanishing baryon-density, giving Tc = T/ ~ 167 MeV for 7s = 1 and 
175 MeV for 7^ = 0.5. 



We have thus specified Tf 
through nucleon percolation. 



Tc at /i = through vacuum percolation and yU/ at T = 



To define a freeze-out curve in the whole T—fi plane, we have to combine the resonance 
gas aspects at low baryon density with the baryon repulsion nature at high baryon density. 
If we move from /x = to finite /x, the resulting medium is no longer an ideal gas, since 
now baryon interactions are present, which are not accountable in terms of resonances. To 
reach freeze-out, the system has to expand and cool off enough to stop these non-resonant 
baryon contributions. It is thus clear that for /x 7^ 0, freeze-out will occur for Tf < Tc, 
below deconfinement. In particular, as is seen from Fig. |21 when we reach fi = fif, the 
interacting hadronic medium formed at the confinement point still contains mesons; but 
these have disappeared when the medium has cooled off enough to freeze out all dynamical 
baryon repulsion, leaving cold nuclear matter of standard density. 

Given the overall hadron density n(T, /i) and the net baryon number density nB{T, fi) 
through Eq'ns. (fTTj) and (US)), respectively, we can specify a freeze-out temperature show- 
ing the behaviour just outlined. We assume that the sector of vanishing baryon density 
freezes out according to the resonance gas approximation and by vacuum percolation, 
that of finite baryon number according to baryon percolation, to obtain 



n(T, /i) 



1.24 



1 



nB{T,fi) 
n{T, /i) 



+ 



0.34 



niT^jj) 



(16) 



as the defining equation for the freeze-out curve. It is clear that when /i = 0, we recover 
condition (0), while in the limit of a cold nucleon gas, with n/ns = 1, we get back Eq. 

Moreover, it determines freeze-out fully in terms of geometric clustering based on the 
intrinsic hadronic scale. The model contains no adjustable parameters, with the exception 
of a possible 7^ < 1. In principle, 7^ could depend on the initial collider energy or on fi. 
However, variations of 7^ have only little effect on the resulting curve, as already noted 
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for the limiting case of /i = 0, and this effect decreases further for higher /i. Therefore we 
will perform our calculations for fixed 7^ = 0.5. 

We now turn to a comparison of our results to the relevant data. Starting with data 
for freeze-out temperatures at /i = 0, we show in Fig. El a recent compilation of data from 
elementary hadron collisions 4J; it is seen that the range T ~ 170 ±5 MeV obtained from 
the vacuum percolation condition (jH) for 0.5 < 7s < 1 is in good agreement with the data. 
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Figure 3: Freeze-out temperature in elementary hadron collisions, from 0. 



200 

180 
160 
140 
^ 120 

CD 

^ 100 

I 80 
60 
40 
20 



^, MeV 

Figure 4: Freeze-out in heavy ion collisions, from jTTj. 

Coming to heavy ion collisions, we show in Fig. 0] results from the four different fa- 
cilities. Also here the agreement is seen to be very good, with clear deviations only for 
the lowest energy data from SIS. While we have no explanation for this, it should be 
noted that at such low energies, secondary hadron production is a rather rare process, 
making a freeze-out determination from species abundances more difficult. In particular, 
the number of nucleons is now fixed by the mass of the colliding nuclei, not by any ther- 
mal parameter. It is rather surprising that the data fall below the expected curve. Our 
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assumption of an ideal nucleon gas at high baryon density could well be an oversimplifi- 
cation; but we saw above that further dynamical repulsion, beyond the effect due to the 
Fermi statistics of the nucleons, in fact leads to a freeze-out at higher /i. One might thus 
wonder whether threshold effects at SIS energies could lead to an effective reduction of 
the freeze-out temperature. 

Finally, we compare our freeze-out condition with that obtained from < E > / < 
N >= 1 GeV [7j. In Fig. we show both forms together with the data. While the two 
models never differ by more than 10%, the form of 0, in contrast to ours, provides good 
agreement also for the SIS point. On the other hand, < E > / < N >= 1 GeV leads to 
a rather sudden change of behaviour for fi > 0.85 GeV, which puts nuclear matter below 
the freeze-out point in /x. 
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Figure 5: Freeze-out from percolation and from (E) / (N) = 1 GeV. 
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